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Abstract 

A search for the neutron-unbound nucleus ^^C was performed via the single proton removal reaction 
^ ■ from a beam of ^^N at 68 MeV/u. Neutrons were detected with the Modular Neutron Array (MoNA) 

^ \ in coincidence with ^^C fragments. No evidence for a low- lying state was found, and the reconstructed 

•— H . ^^C+n decay energy spectrum could be described with an 5- wave line shape with a scattering length limit 

^1 of lag I < 2.8 fm, consistent with shell model predictions. A comparison with a renormalized zero-range 

jH . three-body model suggests that ^^C is bound by less than 70 keV. 
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l/^ 1. Introduction 



Since the discovery of the large matter radius of ^^Li pb], neutron halos have been a topic of intense 
^^ \ study near the neutron drip line. The halo structure results from one or more valence nucleons being loosely 

^T) ■ bound which, combined with the short range of the nuclear force, allows them to have a large probability 

of being found at distances much greater than the normal nuclear radius [2|. For two- neutron halos the 
two-body subsystems are typically unbound [3| and knowledge of the basic properties of these subsystems 
. . , is critical for the understanding of the halo nuclei [3|, ^ . 

rS ■ One two- neutron halo candidate is ^^C, consisting of twice the number of protons and neutrons of ^^Li. 

It has attracted significant experimental [5|, |6| and theoretical [7|, |8|, |9|, [lOl, luj attention recently. ^^C was 
first observed to be bound in 1986 [12|. It is a Borromean nucleus because the two-body sub-system ^^C 
had been shown to be unbound jl^, contradicting an earlier measurement which had claimed that ^^C was 
bound [M]. It took nearly 25 years before first properties of ^^C were measured. A large matter radius 
was extracted from the measured reaction cross section, suggesting that ^^C exhibits a two-neutron halo [5[. 
Using a simple model relating the measured radius and the two- neutron separation energy S2ni the authors 
deduced a strong s-wave configuration in ^^C. Subsequent momentum distribution measurements following 
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Figure 1: Secondary beam composition with -^-^N beam of interest indicated along with primary contaminants. The left panel 
shows the identification by energy-loss and time-of-fiight. The shaded region in the right panel shows the applied gate on the 
time-of-fiight parameter to select ^^N. 



neutron removal reactions supported this suggestion of a significant i^s'^,^ valence neutron configuration 
in ^^C [6]. Very rece ntly the mass excess of ^^C was measured to be 53.64(38) MeV corresponding to 



S2n = —140(460) keV [15|. Given the fact that ^^C is actually bound this value probably should be better 
quoted as S2n = OIq^^ keV. Theoretically, an upper limit of 220 keV for S2n was calculated assuming a 
dominating 5- wave contribution [9|. 

These observations strongly support the presence of a two- neutron halo in ^^C. This typically implies that 
the "two-body subsystems must be either very weakly bound or low- lying resonances, or virtual states must 
be present, in order to support a halo state" [3]. As mentioned earlier, ^^C is unbound but no spectroscopic 
information is presently known. The ground state of ^^C is expected to be a 1/2+ state with a large i^sl^2 
single particle configuration although the possibility of a degeneracy or even level inversion with the i^d^/2 
has been suggested [l6| . In the present paper we present a search for a low lying resonance/virtual state in 
^^C. We attempted to populate ^^C with one-proton removal reactions from a secondary beam of ^^N and 
extracted the decay-energy spectrum via invariant mass spectroscopy by measuring neutrons in coincidence 
with 20c. 



2. Experimental Setup 

The experiment was performed at the National Superconducting Cyclotron Laboratory (NSCL) at Michi- 
gan State University. A 90 pnA ^^Ca primary beam at 140 MeV/u impinged on a 2068 mg/cm^ ^Be pro- 
duction target, and isotopic separation of a 68 MeV/u ^^N secondary beam was achieved using the A1900 
fragment separator [17] with a 1057 mg/cm^ Al achromatic wedge degrader placed at the dispersive image. 
All data were taken with A1900 momentum slits set at 2.5% acceptance, which resulted in a ^^N particle 
rate of 37/s with a purity of 32%. The beam composition is shown in Figure [H The two heavy ion contam- 
inants were ^^F at 6% and ^^C at 2.8%; light ions comprised the rest of the contamination. The energy loss 
(AE) versus time of flight plot shown in the left panel demonstrates the clean separation of ^^N from the 
contaminates as a function of time-of- flight. The shaded area in the right panel shows the applied gate in 
time-of-flight to select ^^N. 

Figure [2] shows a diagram of the experimental setup. After exiting the A1900, the secondary beam passed 
through two position-sensitive cathode readout drift chambers (CRDC-1 and CRDC-2). Downstream of the 
CRDCs, a quadrupole triplet magnet focused the beam onto the 481 mg/cm^ ^Be reaction target. A 
0.254 mm plastic scintillator (SCI-1) was placed immediately upstream of the reaction target to determine 
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Figure 2: Diagram of experimental setup. 

beam, charged fragment, and neutron time of flight. The flight path length from the A1900 to the target 
position was 11.6 m. Neutron unbound isotopes produced in the reaction target immediately decayed into 
charged fragments and one or more neutrons. A large-gap superconducting dipole magnet [18] bent the 
charged fragments away from the beam axis, and the neutrons were detected near zero degrees by the 
Modular Neutron Array (MoNA) [Toj . 

The magnet has a maximum rigidity of 4 Tm, a vertical gap of 14 cm through which neutrons pass, and a 
bending angle of 43°. The magnet was set to a rigidity of 3.8 Tm to center the ^^C fragments on the charged 
particle detectors. Downstream of the magnet, two 30 cm square CRDCs separated by 1.82 m provided 
fragment trajectory information (CRDC-3 and CRDC-4). An ionization chamber and 4.5 mm thick plastic 
scintillator (SCI-2) provided energy loss information for element separation, while a 150 mm thick plastic 
scintillator (SCI-3) provided a total residual kinetic energy measurement. 

MoNA [19] consists of 144 10 x 10 x 200 cm^ plastic scintillator bars with photomultiplier tubes (PMTs) 
attached to each end. The modules were arranged in walls that were 16 modules tall and centered on the 
beam axis. Walls of 2 by 16 modules each were positioned with their front faces at 5.90 m, 6.93 m, and 
7.95 m from the reaction target. A block of three walls was placed at 8.65 m. Neutron time of flight was 
measured from the mean time of the two PMT signals of the detector module that detected an interaction, 
while position across the bar was measured by the time difference between the signals. 

3. Data Analysis 

A one-proton removal reaction was used to populate neutron unbound ^^C from the ^^N secondary beam. 
Charged fragments were separated by element using energy loss information from the ionization chamber and 
the time of flight between SCI-l and SCI-2 as shown in Figure [3l Finally, isotopes of a given element were 
separated by correcting the time of flight between SCI-l and SCI-2 for correlations with both dispersive and 
nondispersive angle and position measured by CRDC-3 and CRDC-4. This separation is shown in Figure HI 
with the grey region representing the selected ^^C fragments of interest. The solid line indicates the result 
of a fit including three convoluted Gaussian functions to estimate cross contamination between isotopes, 
while the dashed lines indicate each isotope's contribution to the total. With the ^^C selection gate shown, 
the contamination level is 2% and 20% of the ^^C fragments are lost. 

Reconstruction of the decay energy of ^^C was performed using the invariant mass method. In this 
method, the decay energy E^ can be expressed as 
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where Ef and En are energies, pf and Pn are momentum vectors, and rrif and m^ are masses for the 
charged ^^C fragment and neutron respectively. 
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Figure 3: Element separation via energy loss in the ionization chamber vs. time of flight. 



Energies and angles for the charged fragments at the target were constructed from the post-magnet 
trajectory information using a transformation matrix generated by COSY INFINITY [20| based on measured 
magnetic field maps. The reconstruction was improved by including the measured dispersive position at the 
target in the transformation matrix as described in Ref. |21|. The dominant uncertainties for this process 
were position and angular resolution in CRDC-3 and 4, with values apos = 1.3 mm and aang = 1 mrad 
respectively as determined by data taken with a tungsten mask placed in front of the detectors. 

Neutron energies and angles were determined from the position and time of fiight of the neutron inter- 
actions in MoNA. The time of fiight resolution had atof = 0.3 ns, while position resolution along a bar was 
characterized by shadow bar data and simulations as a sum of Laplacians 22|, |23| 
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where / = 0.534, ai = 16.2 cm, and (72 = 2.33 cm. Discretization of the detector bars in the y and z 
directions resulted in uniform uncertainty of ±5 cm for those components. 

Monte Carlo simulations of the full experimental setup, reaction process, and decay characteristics of 
^^C were performed to account for the effects of each detector's geometrical acceptance and resolution on 
the reconstructed decay energy. The reaction in the target was treated by a Goldhaber model [2J with a 
friction term [25| . These simulations were validated and all free parameters were constrained by comparison 
to experimental spectra except for those related to the neutron decay. Figure [5] shows the comparison 
between simulation and data for fragment position in CRDC-3 and reconstructed fragment kinetic energy. 
Simulated data sets were generated to populate grid points of the decay parameter phase space in question; 
these were analyzed using the same software as the experimental data and the two data sets were directly 
compared. 



4. Results 

Figure [6] shows the measured decay energy for ^^C. The data (black squares) are distributed over a broad 
energy range and do not exhibit any obvious resonance. A sharp resonance was not expected because the 
i = 2 states were not expected to be populated in the proton removal reaction from the ^^N ground state, 
which is accepted to possess a J^ of 0~ 26|, |27|. The calculated spectroscopic factors for the £ = 2 5/2+ and 
3/2+ excited states in ^^C are and 0.05, respectively. In contrast, the spectroscopic factor for populating 
the ^ = 1/2+ state is 0.75. 
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Figure 4: Particle identification for neutron coincidence data showing 18-20 q with dashed fines demonstrating the relative 
contributions of the individual isotopes. The shaded region shows the selection of -^^C, with an estimated 2% contamination 
from -"^^C and 20% losses from the gating procedure. 

Therefore, the decay energy spectrura was fit assuraing a pure s-wave decay, with the line shape for the 
decay calculated from the analytic approximation from Eq. (30) of Ref. |3l| : 
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where the parameterization S = ask is used for the phase shift, 7 = ^/—2mei/fi is the decay length of 
the initial ^^N ground state, e^ is the binding energy of ^^N, k = s^2mef /% is the final momentum of the 
neutron in the continuum, and e/ is the decay energy of the neutron. 

The scattering length a^ was allowed to freely vary in the fitting process from to —100 fm. The 
fitting process used a minimum x^ technique and favored a la limit of \as\ < 2.8 fm, with the best fit 
being as = —0.05 fm (shown in Figure [6] as the dash-dotted and dotted lines respectively). The figure also 
demonstrates the sensitivity of the present set-up to any potential low energy virtual state by comparing an 
5- wave decay with a^ = — 15 fm (dashed line) to the data and short scattering length curves. It is clear that 
any low lying 5-wave state would have been apparent in the decay energy spectrum. The small magnitude 
of the extracted scattering length of a^ = —0.05 fm is essentially equivalent to no interaction at all, thus 
quoting a resonance energy is not meaningful |32J . 

It should be mentioned that we cannot rule out the possibility of decay to the 2+ state in ^°C as the 
setup did not measure coincident 7-rays. However, it is not likely that this decay mode be observed. The 
spectroscopic factor for decay to the 0+ ground state of ^^C was calculated to be 0.38, while that for the 
2+ first excited state was 1.2. Since the spectroscopic factors are relatively close in magnitude, the lack 
of momentum barrier for ^ = neutron decay to the 0^ will cause it to dominate the i = 2 transition to 
the 2+. Furthermore, the £ = 2 transition to the 2+ would appear as a well-defined resonance in the decay 
energy spectrum, and no such resonant peak is observed. 



5. Discussion 

The results of different shell model calculations for ^^C and ^^C are shown in Figure [71 together with 
the experimental level scheme for ^^C and the recently measured S2n for ^^C. We performed calculations 
with NuShellX@MSU [33] in a truncated s —p — sd — pf model space with a modified WBP [3J interaction 
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Figure 5: Comparison of data (squares) with simulation (lines) for (a) ^^C fragment position and (b) 
energy after the reaction target. 
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labeled WBP*. In this interaction the neutron sd two body matrix elements (TBME) were reduced to 
0.75 of their original value in order to reproduce a number of observables in neutron rich carbon isotopes 
[29[. Comparison was then made to the unmodified WBP interaction as well as the WBT and WBT* 
interactions, where WBT* incorporates the same TBME reduction as WBP*. These interactions have 
been used extensively for calculations in this neutron rich region with Z<8 [l6|, |26|, |29|, |35|, l36|. For all 
interactions, the ^^C 1/2+ s-state is expected to be unbound by more than 1.4 MeV, which is consistent 
with our measurement as no low lying states are predicted. However, as seen in the figure, all interactions 
also incorrectly calculate ^^C to be unbound. 

Two recent shell model calculations with different interactions do not have the same problem. Coraggio 
et al. derived the single-particle energies and the residual two-body interaction of the effective shell-model 
Hamiltonian from the realistic chiral NN potential N^LOW [30|. The calculations predict ^^C to be unbound 
by 1.6 MeV and ^^C to be bound by 601 keV which is even more bound than the recent measurement of 
0^0^^ keV by Gaudefroy et al [l5[. The calculations by Yuan et al use a newly constructed shell- 
model Hamiltonian developed from a monopole-based universal interaction {Vmu) PjJ- The results of the 
calculations are consistent with the two-neutron separation energy of ^^C and predict ^^C to be slightly less 
unbound (~1.3 MeV, extracted from Figure 6 of Ref. [ll|) as compared to the other calculations. 

All models predict ^^C to be unbound by significantly more than 1 MeV. This is consistent with the 
non-observation of any resonance structure in the present data. However, it is somewhat unexpected as 
it represents the first case where the two-body subsystem of a Borromean two-neutron halo nucleus does 
not have a low-lying virtual state or resonance below 1 MeV. This might be another indication that the 
two-neutron separation energy of ^^C is extremely low and very close to the binding threshold. 

The relationship between the two-neutron separation energy of a halo nucleus and the virtual state or 
resonance of the unbound two-body subsystem has been demonstrated before for ^^Li [4]. Higher resonance 
energies or smaller absolute values of the scattering length of the virtual state correspond to smaller two- 
neutron separation energies. Yamashita et al. derived this relationship for the ^^C-^^C system within a 






50 


1 

0.8 



5 0.6 

N 



gO.4 



0.2 




Gs = —15 fm 
Gs = —2.8 fm 
a, = -0.05 fm 
data 



vH^ 



0.5 1 1.5 2 2.5 

Reconstructed Decay Energy (MeV) 



Figure 6: Reconstructed decay energy spectrum (including experimental acceptance and resolution) for ^-"^C where squares 
represent data points. The dotted, dash-dotted, and dashed lines represent s-wave curves of as = -0.05 fm, -2.8 fm, and -15 fm 
respectively. The inset shows the RMS decay energy resolution of the setup as a function of decay energy, and the top panel 
shows geometric efficiency of the setup as a function of decay energy. 



renormalized zero-range three-body model assuming a pure s-wave valence neutron configuration 7, |8[. 
They establish limits on S2n of '^120 keV and ~70 keV based on the measured matter radius of ^^C [5| for 
virtual state energies of keV and 100 keV, respectively. 

For states near threshold the energy of a virtual state can be related to the scattering length with 
E c^ h /2/ia^, where /i is the reduced mass [37|. This relationship places a 100 keV virtual state in '^^C at 
approximately a^ = — 15 fm. As shown in Figure [6] our data are not consistent with such a low- lying state, 
limiting the two- neutron separation energy of ^^C to less than 70 keV. 



6. Conclusions 

In summary, we have searched for the neutron-unbound nucleus ^^C via one-proton removal from a 
^^N beam. The reconstructed ^^C+n decay energy spectrum does not contain any evidence for a low- 
lying state, which agrees with shell model calculations that predict this state to be neutron unbound by 
1.5 MeV. From this non-observation an upper limit on the two- neutron separation energy of ^^C can be 
placed at 5'2n < 70 keV based on calculations from a renormalized zero-range three-body model [71, |8[. 
This is consistent with the recently measured limit of 5'2n < 320 keV [15| as well as calculated limits of 
6'2n < 220 keV ^ and 6'2n < 50 keV ^. 

In the future, the search for ^^C should concentrate on populating the ^5/2 state. Although the Si/2 is 
still predicted to be the ground state it cannot be observed at the calculated decay energy. The ^5/2 state 
can be populated with a one- neutron transfer reaction with a ^^C beam on a deuteron target. The excitation 
energy is predicted to be above the first excited 2+ state of ^^C and thus 7-rays have to be detected in order 
to determine if any observed decay populates the ground state or the first excited state of ^^C. 
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Figure 7: Experimental levels of ^"-"C and theoretical 1/2+ and 0"^ ground states for "^-^C and ^^C, respectively. The data for 
^°C are from [281, [29l and the measured two-neutron separation energy for ^^C from [l^. The N'^LOW and VMU calculations 
are from [SQl and [ll|], respectively. 
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